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Abstract 

General solutions of nonlinear ordinary differential equations (ODEs) are in general difficult to 
find although powerful integrability techniques exist in the literature for this purpose. It has been 
shown that in some scalar cases particular solutions may be found with little effort if it is possible 
to factorize the equation in terms of first order differential operators. In our present study we use 
this factorization technique to address the problem of finding solutions of a system of general two- 
coupled Lienard type nonlinear differential equations. We describe a generic algorithm to identify 
specific classes of Lienard type systems for which solutions may be found. We demonstrate this 
method by identifying a class of two-coupled equations for which the particular solution can be 
found by solving a Bernoulli equation. This class of equations include coupled generalization of the 
modified Emden equation. We further deduce the general solution of a class of coupled ordinary 
differential equations using the factorization procedure discussed in this manuscript. 

PACS numbers: 
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I. INTRODUCTION 



The modified Emden equation (MEE) 

u + auii + (3u 3 = 0, (1) 



where overdot denotes differentiation with respect to time, has been studied by Painleve 



more than a century ago. He and his coworkers found that ([I]) is integrable for two specific 
parametric choices, namely (i) f3 = —a 2 and (ii) /3 = Recently Chandrasekar et al ^ 
have shown that this equation is integrable for any choice of arbitrary parameters a and 
(3. The same authors have also underlined its immense applicability through its intimate 
connections with other well-known systems like the force-free Duffing-type oscillator and the 
two-dimensional Lotka-Volterra system. It also has substantial importance in physics itself 
as it rears its head in a plethora of varied contexts such as the equilibrium configurations 
of a spherical gas cloud acting under the mutual attraction of its molecules and subject to 
the laws of thermodynamics (3-6], and in the modeling of fusion of pellets Q. It is also the 
or the spherically symmetric expansion or collapse of a relativistically 



governing equation 
gravitating mass 



-llOj. Equation ([I]) can also be seen as a one-dimensional analog of the 



boson 'gauge-theory' equations introduced by Yang and Mills [111 . Il2| . This equation also 
comes up in a variety of mathematical problems such as univalued functions defined by 
second-order differential equations [l3| and the Riccati equation fill ]. 

Equation ([I]) is a particular case of a more general family of equations of the Lienard 

type, 

u + f(u)u + g(u) = 0. (2) 

The Lienard family of equations is in general not integrable, and only for some specific choices 
of the arbitrary functions f(u) and g(u ) complete solutions can be found. Its applications and 



properties have been discussed in |14j . Il5j . However, in general, the task of finding general 
solutions of such nonlinear second order systems requires the use of complicated, albeit 
powerful, procedures like the Painleve analysis, modified Prelle-Singer procedure, use of Lax 
oairs and the associated inverse scattering transforms. Recently Rosu and Cornejo-Perez 



16 



16| have suggested that at least in the case of some polynomial nonlinearities particular 



solutions may be found rather simply by an elegant method of factorizing the system (J21) into 
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two first order differential operators. Berkovich 18j has explored the factorization technique 
for solving scalar ODEs in substantial detail. The factorizing procedure is also used to 
obtain the travelling wave solutions of certain nonlinear partial differential equations such as 
the Korteweg-de Vries-Burgers equation, Camassa-Holm equation, Kolmogorov-Petrovskii- 



Piskunov equation 
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19| . However this approach has not yet been applied to the study 



of coupled systems. Our motivation in this present study is to extend this approach to the 
study of coupled systems and identify which systems allow solutions to be found by such 
means. 

In our present paper we study a system of two-coupled Lienard type equations, 



u + v)u + f 2 (u, v)i) + #i(it, v) = 0, 
v + f 3 (u, v)u + / 4 (w, v)v + g 2 (u, v) = 0, 



(3a) 
(3b) 



and investigate whether specific forms of this generic equation allow solutions to be found by 
the factorization technique. Here we show that this can be done by factorizing the system 
OH]) in terms of first-order differential operators by generalizing the procedure of Reyes and 
Rosu 20] for scalar second-order ODEs. We find that for the following general class of 
equations finding a solution is as simple as solving a Bernoulli equation in one variable: 

u - b 2 T]ipi + &i£-02 + (hfa _ b 2 ipi)h + (hip 2 - b 2 ipi)h + S[b 2 (r] - di)(h + di)ipi 

-b 1 (£-d2)(h + d 2 )M = 0> ( 4a ) 
v + a 2 rjip l - - (ai^2 - a 2 tp 1 )h-(a 1 'tl) 2 - a 2 ipi)h - S[a 2 (r] - di){h + di)^ 

- ai(£ - d 2 ) (h + d 2 )ip 2 ] = 0, (4b) 

where ^1,2 — a>i,2"u + b\^ 2 v + Ci^ 2 , 77, £ are arbitrary functions of u and v, h is an arbitrary 
homogeneous function in ipi and ip 2 , 5 = a±b 2 — a 2 b± and dx >2 , ai j2 , bi >2 , ci >2 are all constants. 
For suitable choice of arbitrary functions and parameters, namely rj — h, £ — h, tpi — 
u, ij) 2 = v and h = —(kiu + k 2 v ), where ki and k 2 are constants, this reduces to the form 



u + 2{kiu + k 2 v)u + (kiii + k 2 v)u + (A^-u + k 2 v) 2 u + ufu = 0, 
v + 2{kiu + k 2 v)v + {k\ii + k 2 v)v + [kiu + k 2 v) 2 v + u 2 v = 0, 



(5a) 
(5b) 



where uj\ 2 = d\ 2 . This system has been studied in the literature 



2l| and shown to be 



completely integrable by Chandrasekar et al [22]. Also (jSJ) is a two-coupled version of the 
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special modified Emden type equation with additional linear forcing, 

u + 3kuu + k 2 u 3 + uj 2 u = 0. (6) 



23]. 



This equation exhibits the isochronous property 

The plan of the paper is as follows. In Section [Til we investigate the factorization technique 
for the two-coupled system ([3]) in terms of first order differential operators. In section II III 
we introduce a generic algorithm to identify classes of coupled Lienard type equations (j3J) 
whose solutions can be found by simple methods and chalk out the procedure to identify 
a class whose solution may be found in terms of the Bernoulli equation. In section IIVI 
we discuss the procedure to deduce the general solution of a class of coupled ODEs using 
the factorization procedure. We illustrate the procedure by considering a specific equation. 
Finally in section [V] we conclude by summarizing our main results. In appendix A we prove 
that the functions ^i >2 in Eqs. (jlj) should be linear in u and v. In appendix B, we discuss 
the factorization procedure for the equations fl23|) in which hi is a function of ipx alone and 
h 2 is a function of if) 2 only. In appendix C, we show that the procedure to obtain the general 
solution for coupled ODEs discussed in this paper can also be straightforwardly used to 
obtain the general solutions of the scalar equations discussed in Ref. [161 ] . 

II. FACTORIZATION OF THE GENERAL CASE 

If the coupled Lienard system f)3aH3 bp can be factorized in the form: 



[D-<f> 1 (u,v)][D-<l> 2 (u,v)]ik(u,v) = 0, (7a) 
[D-fa(u,v))[D-fa(u,v)]iP 2 {u,v) = 0, (7b) 

where D — 4, i[)\ j2 are arbitrary functions of u and v, and fa are functions of u and v to 
be determined, then the problem of finding the general solution of ([3]) may be addressed by 
simultaneously solving the following set of first order differential equations: 

[D-fa(u,v)]P 1 (u,v) = 0, (8a) 

[D-fa{u,v)}P 2 {u,v) = 0, (8b) 

[D — fa(u, v)]ipx(u, v) = Pi(u, v), (8c) 

[D — fa(u,v)]ifj 2 {u,v) — P 2 (u,v). (8d) 
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Also, as pointed out by Rosu and Cornejo-Perez [l6| for the scalar case, the problem of 



finding particular solutions to Eq. (J3J) can be addressed by solving the reduced equations, 

[D — (j>2(u, v)]ipi(u, v) = 0, (9a) 
[D-Mu,v)}Mu,v) = 0, (9b) 

which may in some cases be relatively simple. 

We are thus motivated to explore the question: For what forms of fi and gi in (TJ|) can the 
system of coupled equations be explicitly factorized in the form |?p ?. To answer this question 
we proceed as follows. 

By expanding (j7|) and manipulating the resulting equations appropriately, we get the 
equivalent of (EJ) as 



+ fiu + f 2 v + §i = 0, (10a) 

V+\[(lpluu^2u ~ ^2uu1plu)u 2 + {lpluv^2u ~ ?p2uv4>lu)uV + (V>lw^2u ~ ^2vv^\v)v 2 \ 


+ hu + Uv + g 2 = 0, (10b) 
where the functions 

fl = ^[(03 + (p4)^2u^lv - (01 + fc^lufav ~ 4>2u^l^2v + fau^ly], (Ha) 

h = 7 [(03 + 04 - 01 - 02)^11,^21, - 02iA^2i, + 4>ivi>2^1v], (Hb) 


U = g[<ha^l1p2u - (03 + 04 - 01 - <p2)^lu^2u ~ fau^lu], (He) 

U = 7 [(01 + 2 )^2«^li, - (03 + (Pd^lufav + 4>2v^l^2u - 04i,^2^lJ, (Ud) 




gi = - [0102^^1 



92 = 7 [0304^1^^2 

and the quantity 

& = i>lui>2v - 1p2u4>lv ^ o. 



0304^1^^2] , (12a) 
<k<fahM , (12b) 

(12c) 
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Comparing Eq. ([3} with Eq. (jTOj) . we find ip\, 2u u = i>i,2uv = ^i,2vv — 0. This implies that 
one can assume 

^1,2 = di )2 w + bi^v + ci i2 , (13) 

without loss of generality. Here <2j, bi, and , c« are arbitrary constants. Using ( 1T31) in (I12cf) . we 
infer that 5 = a\b 2 — fl 2^i — constant. Further, defining the variables = <5/j, z = 1,2,3,4 
and gj = 5c/j, j = 1,2, we obtain the relations 

+ hg 2 = 0102-01) a 2#l + hg 2 = 030402- (14) 

From the above two equations in (TI4"]) one can solve for 0! and 03 as 



0i 



Q2^1 + hg 2 

i>2 



(15) 









2 


d 
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dv 


2 
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2 


d 
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2 



Using (fT5"|) to eliminate 0i and 3 in (ITT]) one can, after some manipulations, get a set of 
PDEs which determine 02 and 04: 

-ai(ai#i + big 2 ) - (a x /i + bxh)<Mu (16) 

+ M2) - (ai/ 2 + hf^faipi, (17) 

-a 2 (a 2 ^i + 6 2 # 2 ) - (a 2 /i + b 2 h)^u (18) 

-b 2 {a 2 g x + 62^2) - (a 2 / 2 + hf^foipi. (19) 
Compatibility of f[T5"j) with f|T7|) and of (fT8"j) with f[T9"j) requires that 

— [ai(ai#i + 6^2) + (ai/i + 61/3)^2^] = ^ [6i(aifl'i + &ifls) + (ai/2 + fei/O^i] > 

— [a 2 (a 2 #i + 62^2) + (a 2 /i + 62/3) 02^] = 7^ [&2( a 20i + hg 2 ) + (^2/2 + &2/4j020i] , 
which implies 

^1/2^ - ^1/1 + ( a i^2 - 61/1 + [/au - /u#i)02 + fei^iu - ai^iti = 0, (20a) 
^2/4^ - ^2/3^ + (02/4 - ^2/3 + [/4« - /3^]02)04 + &2#2« - a 2 g 2v = 0, (20b) 



where 



91 = am + h92, 
fi = ai/i + 61/3, 

A = 02/1 + 62/3, 



#2 = a 2 ^i + b 2 g 2l 
h = ai/ 2 + &1/4, 
/4 = 02/2 + 62/4- 



(21a) 
(21b) 
(21c) 



Note that in Eq. (|2T|) . the various functions fas, i = 1,2,3,4 and ^s, i = 1,2, are just the 
coefficients occurring in the coupled ODE (j3J). 

We thus have two first order linear PDEs (j20aj) and (120b[) for 2 and ^4, respectively, 
which may be solved to find the explicit forms of fa and fa. When this is found one can find 
the explicit forms of fa and fa from (j!5p . which imply that the factorization ([7j) is complete. 

Thus we have answered the question posed at the beginning of this section. The system of 
coupled Lienard type equations ([3]) can be factorized in the form (J7|) if and only if the pair of 
first order PDEs (l20l) can be solved explicitly for fa and 4>4 in terms of u and v. Conversely, 
for every solution fa and 4 for (TSUI) , there exists a factorization (j7j)) corresponding to a 
system of coupled second order nonlinear ODEs of the form fl3]). 

III. OUTLINE OF THE PROCEDURE 

We next turn to the problem of finding solutions by using the factorized formQ. If the 
lower order forms (jSJ) and can be associated with some known system of equations for 
which the solution is known, then the solution of the original coupled second order system (J3]) 
can be found. With this motivation we explore the inverse problem: Given some particular 
forms of (pi, i = 1,2,3,4, satisfying Eqs. / TTjj) and $EW> what is the subclass of two- coupled 
Lienard type systems / TTOj) that can be factorized to the form |2P and (OJ), for some choice of 
arbitrary constants Oj, 6j, Q ? 

We start off by identifying such an interesting choice of fa. 
Rewriting (j5J) we have 

fa = fa fa, fa = fa fa (22) 

We consider the case that the reduced equations are of the form 

fa = (ui + hi(fa,fa))fa, (23a) 
fa = {u 2 + h 2 (fa,fa))fa, (23b) 

where oj\ and u 2 are constants, so that fa = u 1 + hi^fa^fa), fa = u 2 + h 2 (fa,fa). If hi is 
independent of fa and h 2 is independent of fa, then the reduced equations are uncoupled 
and the problem of finding particular solutions reduces to the much simpler problem of 
solving a pair of scalar first order ODEs. An example of this is dealt with in Appendix |B] 
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Otherwise if hi = ah 2 where a is a constant then one integral can be found by eliminating 

aipiip2 - ip2tpi = (awi - u 2 )4>iip2- (24) 



hi from f)23ap and (123b|) as 



Eq. (l24l) can be rearranged to 



a- = aui - u 2 ). 

^1 ^2 



(25) 



Integrating Eq.fl25l) one gets 

^ = C ie (^-^ 2 E C (^ 2l (26) 

where C% is the integration constant, c(t) = de^ 1 '^ 1 . 

Let us now consider the specific case hi = ah 2 = h(ipi,if) 2 ), where h(ipi,if) 2 ) is a polyno- 
mial containing N terms and is of the form 



N 



h(ipi,tp 2 ) = ^klpi^: 



9, 
2 ! 



(27) 



i=l 



where Pi and % are real constants. Using the relation ( 126]) and rewriting the above equation 
we get 



N 



h = Y,h (c(t)^)- Vf. 



(28) 



i=l 



The above polynomial becomes homogeneous for the condition — + = m, where m is a 
constant, 

N 

h = ^^k iC {t) V -i. 

i=l 

Next, substituting f fl3|) in ( ITUl) and simplifying we get 
1 



(29) 



u + 



1 



1 



1 

'"5 



(a 2 6i0 3 + w 2 + /i) - aib 2 ((j)i + u x + h) + (b x ip2 - b 2 ipi)h u 
bib 2 ((p 3 - 4>i + uj 2 - ui) + (6iV»2 - i^i)^. 

&2<Mk>l + /i)^l - &1 3 (W 2 + /i)^2 =0, 

oi&2(03 + u 2 + h) - a 2 bi((pi +ui + h) + {a^ 2 - a 2 i/ji)h 
aia 2 (0 3 - 0i + w 2 - wi) + (aii[> 2 - a 2 4>i)h u 
-- a 2 gi(ui + h)i>i - aig 2 (u 2 + h)i> 2 = 0. 



(30a) 



(30b) 
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where h is given by f )29|) and 0i and 03 are given by 

01 = —7 , 03 = —7 • 31 

■0102 ^204 

Redefining now the constants and functions as ui = 5d%, u 2 = Sd 2 , 0i = S(rj — 
d\)i 03 = <K£ ~ d 2 ), h = 8h, where 5 = ai&2 — a 2 bi =constant, we find that ( 150]) 
reduces to the coupled Eq. 01]). A particular solution of ([30]) can be obtained by solving the 
following Bernoulli type equation 

^ 2 = u 2 i; 2 + d(t)iJ 2 n+1 , (32) 



where d(i) = YliLi kic(t) « , which is obtained by substituting (129]) in Eq. (I23b|) . Equation 
([32]) admits the following explicit solution 

ifa(u,v) = — , (33) 

[C 2 - mj e muJ2t d(t)dt] m 

where C 2 is an integration constant, and, since ipi = if) 2 c(t), we have 

M% «) = ^ r . (34) 

[C 2 -mJ e muJ2t d(t)dt] m 

Inverting the relation ([T3]) one can obtain u and v as 

u = y(-0i - c x ) - -j(ip 2 - C2), (35a) 
v = y(0> 2 - c 2 ) - y(0i - ci), (35b) 
Using the explicit forms of ifii and -02 as given by Eqs. (I3"3"]) and (13~4"]) . we can then write 

«=T I ~ 7TX - c i I - T f ~ " — X - c 2 T (3G ) 

\ [C 2 - to / e^^d^dt] " 

a 2 f c(t)e W2 * 
[C 2 -mj e^^d^dt] ~™ / 5 I [C 2 - m j e^H^dt] ~~ 

where C% in c(t) (vide Eq. ( 126]) ) and C 2 are arbitrary integration constants. 




ci ,(37) 



A. An example : Case ip\ = u, ip 2 = v 



We now consider a specific equation belonging to the class of coupled ODEs ([301 and 
obtain a particular solution using the above factorization procedure. For illustrative purpose, 
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let us consider the simple case ipi = u and i/j 2 = v for which Eq. (130bj) reduces to the form 



u — uh + [h + u\ + 4>i]u — h v ui)-\-fau(h + uj\) = 0, (38a) 
v - h u vii + vh + [h + u 2 + fajv+favfo + u 2 ) = 0. (38b) 

Choosing the arbitrary functions 

fa — — (kill + k 2 v) - ui, 03 = —(kiu + k 2 v) - lu 2 , 

we find (1551) reduces to the two-coupled version of the modified Emden equation ([5j), studied 



in the literature 



21 



22J. Substituting these forms of ipi, ip 2 and h in (1551) we get the following 



particular solution of (jSJ) 



Cio;ia; 2 e^ f 

(39b) 



where C\ and C 2 are arbitrary constants. One can check that this particular solution can be 
obtained from the general solution (13) given in Ref. |2l| after fixing two of the integration 
constants. 

We are also exploring further how the forms (123]) can be generalized so that more general 
system belonging to the class OH]) can be brought into the above formalism. The results will 
be presented in future. 

IV. METHOD OF CONSTRUCTING GENERAL SOLUTION 

In the previous section we have obtained the particular solutions for a class of coupled 
second order ODEs by factorizing them. In this section we obtain the general solution of 
a subset of these coupled nonlinear ODEs ( flOl) by factorizing them for suitable parametric 
choices. For this purpose we assume fa = <fi + X*> * = 1)2,3,4, where Xii i = 1,2,3,4, 
are arbitrary constants and <f> is some arbitrary function of u and v. Then the equation (JHJ) 
becomes 

[D-(f>(u,v)- X i}Pi(u,v) = 0, (40a) 
[D-<f>(u,v)- X 3\P2(%v) = Qi (40b) 
[D - <f>(u, v) - x 2 ] V>i (u, v) = Pi (u, v) , (40c) 
[D - cf>(u, v) - xdMu, v ) = P ^ v). (40d) 
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This implies that the following relation holds good : 



D 



P, 



Xijf + 1, 3 = 1,2 



(41) 



where Xi — (X2 — Xi) an d X2 = (x<i ~ Xs)- Upon integrating equation (j4"Tj) we find a solution 
of the form 



§ =e &*(j J_ e -*M = i,2 



(42) 



where Jj and /2 are two integration constants. Using ( )42l) . one can rewrite Eqs. f l40cl) and 
fliOdl) as 

Pj _ D[^j(u,v)] . e - ** 



Xj 



, 3 = 1,2 



(43) 



where Xi = X2 and x 2 = Xi- From (14"3"|) we get 

"01 ^2 " 



log 



X2 



1 Xi 



(44) 



Integrating equation (jUJ) we obtain 
where 



1 Xi 

j - ^ e -X2« 

Z X2 



(45) 



a(t) 



/1 - ^e"* 1 * 

1 Xi 

1 X2 



Let us assume that the function <fi(u,v) = hfyi,^) = h(u,v), where hfyi,^) is given in 
(|27|) . Now the equation for -i/^ (see Eq. (143}) ) becomes 



h e"* 1 * 

Xi 



m+l 



(46) 



where 



N 



s(t) = ha(t) q \ m = pi + q u 



fco, fci, . . . , fcjv are arbitrary parameters. 
Equation ( 1461) has the explicit solution 

V>i(«,v) = 



1 XI 



-mxit 



I4 — m J 



1 XI 



e m xit s {t)dt 



(47a) 
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a(t) 



From the relations (fT3l) we find 



m 



h 



h 



M 



XI 



-Xlt 



(47b) 



U = yOl - Ci) - j (ip 2 - Ca), 

u = y O2 - c 2 ) - y(^i - ci), 



(48) 
(49) 



where ipi and ^2 are given by Eq. ( 147| . 5 = 0162 — 0261. Note that the above solution 
contains four arbitrary constants Ii, I 2 , I3, and (with I 2 , I3 appearing in the function 
a(t), see Eq. ( I4"5]) ) so that (I47a[) -( |4TE|) constitute the general solution. 



A. Example 

We illustrate the above procedure by considering the following equation belonging to the 
class of equations given by (150]) . 



u + k 2 uv + (Skiu + 2k 2 v — Xi — X2)u + k 2 uv(k 2 v — xi — X2) - kiu 2 (xi + X2 - 2k 2 v) 



+k t u + X1X2U = 0, 



+k^v 3 + X3X^v = 0, 



which can be factorized as 



(50a) 



v + kivii + (3k 2 v + 2kiu — xs — Xi)v + kiuv{k\u - % 3 - xa) - k 2 v 2 (x3 + Xi — 2kiu) 



(50b) 



[D + (km + k 2 v) - xi] Pi(u, v) = 0, 
[D + (k x u + k 2 v) - x 3 ] P 2 {u, v) = 0, 
[D + (kiu + k 2 v) — X2] u — Pi (u, v) , 
[D + (kiu + k 2 v) - Xa] v — P 2 (u, v). 

We obtain the following relations from the above equations, 



. V \ V 



(51a) 
(51b) 
(51c) 
(51d) 



(52a) 
(52b) 
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where xi — X2 — Xi an d X2 = Xt ~ Xs- Integrating we get 



Pi V xi 



Xit 



^ = e^(l 2 --e-^ 
P2 V X2 



Using these relations in (15 left and ( I51dl) we get 



-xit 



u c 

- + (kiu + k 2 v) -xi = 7 

\ Xj 

- + (k x u + k 2 v) - X2 

v ~ i l2 _± e -X2t 



Xa 



(53) 

(54a) 
(54b) 



Integrating the above system of equations we get 

J_ P -xit-\ 



I 2 - -J-e-x* 



X2 



e (X2-Xl)t+h u = 



(55) 



Substituting for v in (154aj) we get 



d 



h - — e"* 11 
Xi 



w + (£4 + k 2 a{t))u 2 . 



(56) 



We wish to note that the above equation falls under the Riccati equation which upon inte- 
gration leads to the following general solution, 

J_g-Xi< \ pXlt 



u 



h 



XI 



I 4 + / (ix - ^e-xiA e»*(fci + k 2 a{t))dt 
Substituting for it in ([33]) we get 



(57) 



I4 + / (l x - ^e-»*j e »*(fci + k 2 a{t))dt 



(58) 



We note here that the above solution contains four integration constants ii, /21 ^3 and J 4 
and this solution can also be obtained by following the procedure discussed in Ref. [24]. 
We find that the general solution of a restricted class of equations, a subset of the class of 
equations ( )30l) , which can be factorized in the form ( 1401) can be obtained using the above 
procedure. For the case of scalar second order ODEs discussed in Ref. [la], obviously we 
can apply the above procedure straightforwardly. The details are given in Appendix C. It 
should be possible to generalize the above procedure to more general cases than (HQ}, though 
we do not attempt this in this paper. 
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V. CONCLUSION 



In this paper we have identified a system of coupled Lienard type equations which can be 
factorized in terms of first order differential operators. We have shown that a particular so- 
lution of equations belonging to this class of coupled Lienard type equations can be obtained 
by solving a Bernoulli type equation. This generic class of ODEs contains the a coupled 
version of the modified Emden equation which has been recently studied in the literature. 
We have also shown that the general solution of a restricted class of equations can also be 
obtained using this procedure of factorization. Several generalizations of our study can be 
proceeded with, by relaxing the various restrictions mentioned in the present work. These 
are being pursued currently. 

In addition to this, one can extend this procedure to higher order scalar/coupled ODEs 
and obtain their corresponding particular/general solutions for suitable choice of the pa- 
rameters. One can also straightforwardly extend this procedure of factorization to a system 
of N coupled second order ODEs. The factorization of N coupled second order ODEs will 
result in a system of 2N coupled first order ODEs similar to (|40l) . From these 2N first order 
ODEs one can obtain a relation similar to Eq. (|4"T|) with j = 1, 2, . . . , N. Integrating this 
system of equations one can obtain the general solution of the underlying TV coupled second 
order ODEs. 
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Appendix A: Proof of V>i,2u« = ^1,2™ = ^1,2™ = 

Requiring that the coefficients of the higher powers of the derivatives contained in the 



bracketed terms in (llObj) be zero we get 



{^Xuu^2v ~ 4>2uu4>lv)u 2 + (Auv^2v ~ ^2uv^lv)uV + {^i vv 1^2v ~ ^2vv^lv)v 2 = (Al) 
(lpluu?p2u ~ ^2uu^lu)U 2 + {lpluv^2u ~ ^2uv^lu)UV + (tp lvv ^ 2u ~ ^2vv^lu)V 2 = (A2) 

Since u and v are independent the coefficients of u 2 ,uv, v 2 in each equation must individually 
be zero. 



^luui>2v 


= tp2uu1plv, 


i>luu^2u 


= ^2uu'^lu 


1pluv^2v 


= ^2uv^lv, 


i>luv^2u 


= i>2uvi>lu 




= i>2vv^lv, 


^lvv^2u 


= ^2vv^lu 



(A3) 



If i>2uu^2uv^2w ^ then (JM]) implies 

Ipluu Ipluv Iplvv i>lu Iplv 



(A4) 



^2uu fauv i>2vv 1p2u 1p2v 

The last equality in ( 1A4j) implies 5 = ipiufav ~ fau^Piv — which is not permissible in ( llObj) 



as in that case the leading order terms il, v vanish and the resulting equation is first order. 

Thus faun, 4>2uv, 1p2vv = 0. 

If ^2uuiip2uviip2w = but one of if)iuu,4>iuv,4>iw ^ then from Eq (jMJ) vp 2u = ip2v = 
which again implies 5 = which is inadmissible. Thus ifii^uu = 4>i,2uv = ^i,2w = 0. 

Appendix B: Separable reduced equations 

We consider the case where hi and h 2 in (I23all23b|) are of the form 



h = Mi, h 2 = k 2 ^\. (Bl) 
The reduced equations are then a pair of Bernoulli equations in ipi and ip 2 

Ipl = UJilpi + h^ +1 , 

^ 2 = u 2 ip 2 + k 2 i) q 2 +1 . (B2) 
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For the class of systems for which ff22]) is equivalent to fl23aH23bl) with hi and h 2 given by 
( 1B1I) explicit solutions can be found relatively simply by solving (IB2j) . Using the condition 
02 = + kiipi and 4 = u 2 + k 2 i\)\ and (fTS]) in (llObj) . we identify this class of systems to 
be 

1 



u + 



{a 2 b 1 g 2 + u 2 + k 2 (q + l)^f) - 0163(51 + Wi + h(p+ 
1 



+ 



5 



6162(52 - 5i + ^2 - Wi + k 2 (q + 1)^1 - + l)^i) 
1 



+ 



5 



1 

V ~~5 



6251(^1 + hipl)^ - hg 2 (u 2 + k 2 ip^)i> 2 
0162(52 + ^2 + & 2 (o + 1)^1) _ 0261(51 + ux + ki{p + 
01^2(52 - 5i + ^2 - wi + k 2 {q + - + 

0251(^1 + hi(%)if)i - a!g 2 (u 2 + k 2 ip^)ifj 2 



0. 



■it 



0. 



where ^1 and -0 2 are given by (|T3l) and g x and 52 are given by (l3Tj) . 
The solution to flB2j) is found using (I33H34I) to be 



i>2 



hide 



k 2 {C 2 e 



-puiit 



-quilt 



bj 2 



Using (135bj) we get an explicit solution for u and v as 



u 



v = — - 



h 

s 

«2 



fci(Cie 



hide 



-puiit V • 

-puiit P \ 



Ci 



Ci 



61 

5 

Oi 

5 



hide 



-quilt 

UJ 2 



hide 



—quilt 



UJ 2 



C 2 



c 2 



(B3a) 



(B3b) 



where d an d C 2 are integration constants. 



Appendix C: General solution of scalar ODEs 

Using the factorization procedure Cornejo-Perez et al. in Ref. [16j have obtained par- 
ticular solutions of a class of scalar second order ODEs. Now, we show that the procedure 
of obtaining general solution discussed in Sec. [IV] is also applicable to scalar second or- 
der ODEs. We demonstrate this by considering a specific equation discussed in Ref. [l^ 
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and obtain its general solution through factorization for suitable parametric choice. Let us 
consider the following equation 

x + {a - ax p )x + px(l - x p )(x p - b) = 0, (CI) 

which is obtained from the Burgers-Huxley equation [16( using a traveling wave reduction. 
We find that for the parametric choice ft = , " 2 ^ 2 , b = — ^(p + 2) — 1, a = — |(p + 2)(a ± 
(X2 ~~ Xi)j Eq. (ICljl can be factorized in the form 

(D-4>- X i)P = 0, (C2a) 
(D — 4>- X 2)x = P, (C2b) 

where D = (p = jpp^x p — \{xi + X2 + a), Xi and ^2 are arbitrary parameters. From Eq. 
(IC2I) and following the procedure discussed in Sec. [TV] we get 



D (ji) =(X2-Xi)| + 1- (G3) 

Integrating we obtain 

p v (X2-X1) 1 j 

Rewriting Eq. ( 1C2I) and substituting for — from the above equation we get 

D\x] , e (xi-X2)t 

0-X2 = t ; (C5) 



where </> = ( p " 2 ) xP — |(Xi + X2 + a )- Integrating the above Bernoulli type equation we obtain 



j _ e(X!-X2)A e _l( X2 _ Xl+a)t 
X2-X1 / 



X W = 7 7 ; zvp — : vp ( C6 ) 



Similarly one can obtain the general solution for all the other equations discussed in Ref. 
[r| for suitable choice of parameters using the procedure discussed here. 
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